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We consider theoretically two coupled optical waveguides with a varying bar- 
rier height along the waveguides direction. The barrier could be constructed 
by the elongated island with a reduced refractive index (which acts as a po- 
tential barrier), such that in the middle region it splits a waveguide into two 
weakly coupled parts. It is predicted by numerical simulations and analytical 
consideration that the presence of some imperfection of the system parameters 
can cause splitting of injected laser beam and one will observe two intensity 
maximums at the output, while for small imperfections the input and output 
beam intensity distributions will be the same. The switching between two 
regimes could be achieved changing spectral width of the beam or refractive 
index of the island. This nontrivial effect is explained by possibility of tran- 
sitions between the different eigenstates of the system in the region of large 
potential barrier heights. The mentioned effect could be used for all-optical 
readdressing and filtering purposes. 
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1. Introduction 

The system of coupled optical waveguides [1, 2] attracts a steady interest of both theoretical 
and experimental researchers since they allow investigating not only purely optical phenom- 
ena (such as anomalous diffraction and reflection [3]) but also generic problems of nonlinear 
physics (creation, propagation and interaction of nonlinear objects [4, 5, 6, 7, 8, 9], various 
types of instabilities [10, 11], nonlinear bistability [12], etc.). At the moment there exists an 
increasing interest to observe various quantum effects in coupled optical waveguides, in other 
words the people started to think to use those systems as macroscopic quantum laboratory, 
as far as the problem could be described by Schrodinger equation, where the waveguide 
direction plays a role of time. Here we just quote the prediction [13, 14] and observation 
[15, 16] of Landau-Zener tunneling, Bloch oscillations [17, 18], Anderson Localization [19] 
and analogy of Josephson oscillations [20, 21]. However, all of these studies investigate 
the dynamical problems, while to the best of our knowledge no statistical research on the 
mentioned systems has been done yet. In the present paper we try to start filling this gap 
investigating the influence of the value of the imperfection of the waveguide and beam pa- 
rameters on the light propagation through the waveguide with the longitudinal potential 
barrier created by the island with a reduced refractive index (see Fig. 1). 

Particularly, we propose splitting of the injected laser beam due to imperfections in the 
system. Indeed, as far as the waveguides are formed by application of external electromag- 
netic fields (for instance, the elongated island of the reduced refractive index could be formed 
by a control beam perpendicular to the sample plane as in Ref. [22]), in this situation due 
to the photon number fluctuations [23] of control beam the refractive index is not homoge- 
neous in space and time. The time averaged stationary picture observed in the experiment 
could be treated as statistical average over multiple distributions of random inhomogeneous 
component of the refractive index. 

In the recently proposed scenario of beam splitting [22] the refractive index defect is 
placed on the way of the injected signal beam. In other words, moving the defect one can 
split the beam or only change its direction. In this paper we do not propose to change 
the position of refractive index barrier, the splitting effect could be achieved by increasing 
the power of the control beam even in the linear regime of the signal beam. We interpret 
predicted effects in terms of establishment of statistical equilibrium between symmetric and 
antisymmetric eigenmodes of the signal beam. 
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Fig. 1. suggesting experimental setup for observation of self-chaotization process between 
two eigenstates. In the proposed device there is an elongated area with a smaller refractive 
index (n < hq) and spatial length Az splitting the waveguide into two parts in the middle 
region. The light is polarized linearly along transverse y direction. In the upper panel the 
low imperfections of the system are considered and and as a result the intensity distibutions 
across the waveguide (axis x) is the same both at the input and output. Lower graph 
represents the situation for larger imperfections and the picture is different: Instead of a 
single intensity maximum at the output one should observe splitting of the injected beam. 
In the lowest graph it is displayed the input-output intensity distributions in case of large 
imperfections. 

In the figure we present schematics for the observation of the predicted effect. The system 
could be effectively considered as a dynamics in a double well potential with a varying barrier 
height. We propose to inject the beam at the left described by symmetric (with respect to the 
center of the waveguide) ground state wavefunction. Then, as it will be shown below, if the 
parameter of self-chaotization is small (i.e. 6ICAz ^ 27r, where 6JC is imperfection parameter 
defined below and Az is a length of the reduced refractive index island, i.e. potential barrier 
region) one will observe the same intensity distribution at the right, while for larger self- 
chaotization parameters the beam will split into two parts. The same switching effect could 
be achieved keeping imperfection parameter 6JC constant and changing height or length of 
the potential barrier (i.e. parameters of elongated island with a reduced refractive index) or 
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controlling spectral width of the injected signal beam. 

The explanation of the effect is as follows: initially, when the barrier height is zero (at the 
left in Fig. 1), the propagation constant of symmetric ground state wavefunction and the 
one of first excited antisymmetric wavefunction differs from one another and no transition 
occurs, thus without no barrier, or low barrier heights the input will coincide with the output 
in any case. In the region with high barrier heights the propagation constants of symmetric 
and antisymmetric wavefunctions are almost the same, thus even small fluctuations (but 
in case of large /S.z for which 51C/S.Z > 27r) can make transitions between those states and 
finally thermalize the situation. In other words the weights of both states will become the 
same and at the right one will observe this mixed state (with two intensity maximums) 
which will be quite different what we inject at the left. Obviously the latter scenario does 
not take place if imperfections are small, then the transition probability between the ground 
and first excited states will be small and the length of region with nonzero barrier heights 
is not sufficient for self chaotization of the states. Thus one will see the same picture at the 
output as at the input. 

Thus we have two types of control parameters: one is the value of imperfections and 
the second is potential barrier characteristics (height and length). Decreasing imperfection 
parameter one needs more and more barrier height and length in order to observe a splitting 
effect. And contrary, in presence of large imperfections in the system even small defect 
with reduced refractive index could split the injected signal beam. It should be especially 
mentioned that in order not to excite also higher order modes the refractive index barrier 
variation should be smooth. 

2. Theoretical Model and Numerical Simulations 

Let us start the quantitative consideration with the Maxwell's equations in a non-magnetic 
medium without free charges, which could be written for the electric field S{r, t) as (for the 
full details of justification of the given procedure and assumed simplifications we redirect 
reader to the Refs. [4, 24]): 



where the refractive index n = n{r) is allowed to vary in transverse to y plane xz (see Fig. 
1 for the axis directions). Then we may assume a field polarized and homogeneous along 




(1) 
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the direction y, namely the transverse-electric field reads as 

S{f, t) = ey (^(a;, z)e-'^'^'-'''^ + c.c.) , (2) 

where Cy is the unit vector of the ^/-coordinate, the stationary envelope function ^(x, ^) 
depends slowly on the coordinate z, uj is carrier frequency of a beam, k = ootiq/c is a carrier 
wavenumbcr and no is an average linear refractive index. In these definitions the Maxwell's 
equation (1) becomes 




% 0.1 0.2 0.3 0.4 0.5 



Fig. 2. Dependence of the propagation constants and of the first two modes on the 
barrier height v of the double harmonic well potential. As seen, for large v-s the eigenvalues 
E^ and E2 are very close. Insets display the form of linear modes in double well potential 
for zero barrier and large barrier heights. 

Now restricting ourselves to the linear regime, let us assume that 1 — has a form of 

square double well potential (considered before in various physical contexts [25, 26, 27, 28, 
29]) plus some random inhomogeneous part H'{z,x). In numerical simulations we choose 
the points in {x, z) plane giving them the random values and then plotting the surface. We 
run the simulations for the multiple random surfaces averaging afterwards the results. 

Defining ujq as a central frequency of the injected beam and rescaling to the dimensionless 
spatial variables z —>■ zno/2u!oC, x —>■ xuq/ujqC, we arrive to the Schrodinger equation written 
in the form: 

1 9^ 
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where V{x) represents some barrier potential centered at x — 0; v{z) is a height of the 
barrier of double square well potential varying along z direction (see for schematics Fig. 
1) and dimensionless variable x varies from —L to L (across the whole transversal length 
of the waveguide), thus \['(z, — L) = \E'(2;,L) = 0. Obviously, in the realistic situation 
one has a finite well depth resulting in a exponential decay of the wave function at the 
boundaries instead of the vanishing condition, but this will cause insignificant modification 
of the results (see e.g. Refs. [25, 26, 27, 28]) and for simphcity reasoning we will consider 
further infinite barrier heights at both boundaries of the double well potential. Under these 
above conditions the eigenvalues of lowest (symmetric with respect to the point x — 0) and 
first excited (antisymmetric) eigenstatcs for various potential heights v are presented in Fig. 
2. 

Without random perturbative terms and constant v{z) equation (4) is easy to analyze 
and we display dimentionless propagation constants and E'^ for first two linear modes 
versus barrier height v in Fig. 2. As seen, if one starts from the symmetric mode and with a 
zero energy barrier {v = 0) eventually increasing barrier height, the propagation constants 
of symmetric and antisymmetric modes become almost equal and infinitesimal fiuctuations 
can cause the transition between those modes. Such a possibility has been discussed earlier 
[31, 32] comparing the scenarios of emergence of chaos in classical and quantum pendula. It 
has been argued that for large barrier heights the energy eigenstates become degenerate due 
to the presence of fluctuations. Therefore, decreasing the barrier height one will not recover 
the initial mode, instead one will have a mixed state consisting of both modes with equal 
weights. In opposite, if one will not increase barrier height enough to consider the states 
degenerated, the system will recover initial mode. 

Our aim is to consider the similar scenario for optical systems like presented in Fig. 1. 
We will monitor the evaluation of the system for a given random distribution of the potential 
H' in (4) and then average the evaluation results over many different random distributions, 
what is presented in Fig. 3. Thus it is assumed that observed stationary picture is a result 
of averaging over fast time fluctuations of the refractive index distribution in the sample. 
In case of both a) and b) graphs we take the same initial mode for zero barrier height 
= ^l/y/T^ cos {'Kx/2L) and the variation of the potential barrier along z is displayed in 
the inset. In the graph a) the imperfections 5K, are small and initial state recovers, while in 
case of large fluctuations displayed in graph b) one monitors at the output both initial and 



6 



first excited modes with the equal weights. 
3. Analytical Consideration 

For analytical treatment we consider first two real symmetric and antisymmetric 
eigenfimctions with Eigenvalues and £'2 of the operator Ho with constant barrier height 
V. Thus the modes satisfy the following equalities: 

and it is assumed that the following orthonormalization conditions hold: 

L L L 

J {^+fdx = 1; J {^-fdx = 1; J = 0. (6) 

-L -L ~L 

The dependence of Eigenvalues and E2 of these modes on the barrier height v is presented 
in the Fig. 2. In case of zero barrier height the eigenfunctions have trivial form (see the 
inset in the same figure): 

$+ = (i/Vl^ cos {7rx/2L) ; $_ = sin(7ra;/L). (7) 

Presenting now wavefunction ^ as 

z) = Mz)^+{x) + (8) 

and defining the symmetric Tig — Ti-^z, —x) + 'H'{z, x) and antisymmetric Ti'a — T~i'{z, —x) — 
H'{z,x) parts of the random potential surface H' we substitute (8) into the (4). Then 
multiplying on $+ and and then integrating over x we get the following set of equations: 

i^^E,^, + }C^2, i^^E,^l;2 + }C^i (9) 
oz oz 

here £'1,2 = E^ 2+^^1,2 ^-^d we have taken into account the following relations and definitions: 

L L L 

J ^±n',^±dx = 5£i,2, J ^±n'^^±dx = 0, J ^+no^^dx = o, 

-L -L -L 

L L 



J <5+n',^-dx = 0, J <5+n'^^-dx = K. (10) 



7 



The set of equations (9) have been considered long before (see e.g. [25, 26]) and has the 
exact solutions. Particularly, choosing the initial conditions as V'i(O) = 1 ^ind ip2{0) = and 
defining AE = E2 — Ei, ipi^2 = \'4'i,2\ exp (i0i,2) we get the solution as follows: 




Fig. 3. Three dimensional graphs of the intensity distribution describing numerical simula- 
tions on the equation (4). The runs are averaged over random distributions of Ti' = Ti'^ + Ti.'^ 
[see equation (4)]. Particularly, for the graph a) the distribution is taken such that a/ (/C^) ~ 
^/{SW) ~ ^/{SW) ~ 0.0002, while for the b) graph ^/{jO) ^ ^/{6Wj ~ ~ 0.002. 

For the larger fluctuations it is clearly seen two intensity maximums at the output part of 
the graph b). The inset in the graph a) shows the variation of the potential barrier height 
V along the waveguide in both large and small fluctuation cases. 

Now we consider two limiting cases: in the first one corresponding to the small barrier 
heights V, one has strongly nondegenerate situation AE ^ /C, the system remains close to 
the initial state and I?/;!! ~ 1 for any z. In the second case of large barrier heights v for which 
E^ ~ E2 it could be shown that 6E^ — 6E2 ^ as well, thus we have the degenerated 
limit AE and from (11) we get: 

V'l = cos(/Cz)e-*^l^ ^2 = -ism{ICz)e-'^^' (12) 

Thus now the intensity oscillates between the states with a wavenumber K,, which could be 
treated as a chaotization parameter. 
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Assuming now that the large barrier heights region lasts for the distance and then 
follows the zero barrier region, from (8) and (12) one can write the wave intensity distribution 
as: 

|^|2 ^ |e-^-^i(--^-)cos(/CA^)$+ - ie-^-^^("-^")sin(/CA^)$_|2 = 

cos2(/CA^)$^ + sm\lCAz)^^_ - sin(2/CAz) sm{AE{z - A^))$+$_. (13) 

As it was expected, expression for intensity contains interferential term. 

Let us note that, difference between pure and mixed states is analogues to the difference 
appearing in case of interfering coherent and non-coherent light fluxes, respectively. As well 
known, in case of pure states the amplitudes are summed in each point, while in case of 
mixed states one has to sum only the intensities. In other words, the interferential term 
existing in case of pure states will disappear in the mixed states situation by averaging with 
respect to random phase. 

Let us suppose that transition amplitude /C is characterized by mean-root-square fluctua- 
tion 6IC = \/ {K?). We discuss here two limiting cases corresponding to the small 5K,A.z -C vr 
and large 5K,A.z ^ tt, respectively. In the flrst case after averaging with respect to the small 
5K,A.z <^ TT we get: 

(cos(K:^)) ^ 1, (cos2(/C^)) ^ 1, 
(sin(/Cz)) ^ 0, (sin2(/Cz)) ^ 0, 

(*) ^ $+(a;), m^) ^l{x). (14) 

It is easy to see that initial wave function in the "input" is completely revived in the 
"output". In case of large quantity 51CJS.Z ^ tt situation is completely different. In this 
case error incursion 51C up to the values 27r takes place on the distance Az and this leads to 
the lost of information about the systems initial state. After averaging procedure for wave 
function we have (^) = and taking into account expressions 

(cos2(/CA^)) = (sin2(;CA^)) = i, (sin(2;CA;s)) = 0, (15) 

we obtain for intensity the following: 

m^) = \{^\{x)^^^_{x)). (16) 
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So, in this case system's state is mixed state in which symmetric and anti-symmetric 
components are presented with equal probabihties. It is important that, condition of self- 
chaotization 51C/S.Z ^ tt could hold even for small 51C increasing the distance Az. 

4. Conclusions 

Summarizing we have investigated the influence of the system imperfections on the intensity 
distribution along the waveguide direction. The predicted splitting of the beam for large im- 
perfections is explained in terms of the chaotized transitions between the quasi-degenerated 
eigenmodes of the Schrodinger equation in double well potential. Thus one can control the 
splitting process by multiple ways: either changing the imperfection value of the system 
or the spectral width of the injected beam or vary the height and/or length of the barrier 
potential. In presence of large imperfections in the system even small defect with reduced 
refractive index could split the injected signal beam. 
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